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Summary. The linearized Boltzmann equation is considered to describe small spa- 
tial perturbations of the homogeneous cooling state. The corresponding macroscopic 
balance equations for the density, temperature, and flow velocity are derived from 
it as the basis for a hydrodynamic description. Hydrodynamics is defined in terms 
of the spectrum of the generator for the dynamics of the linearized Boltzmann 
equation. The hydrodynamic eigenfunctions and eigenvalues are calculated in the 
long wavelength limit. The results allow identification of the hydrodynamic part 
of the solution to the linearized Boltzmann equation. This contribution is used to 
calculate the fiuxes in the macroscopic balance equations, leading to the Navier- 
Stokes equations and associated transport coefficients. The results agree with those 
obtained earlier by the Chapman- Enskog method. The implications of this analysis 
for application of methods of linear response to granular fluids and derivation of 
Green-Kubo expressions for transport coefficients are discussed. 



1 Introduction 

The derivation of hydrodynamic equations for granular gases and identifica- 
tion of expressions for the associated transport coefficients has been a prob- 
lem of theoretical and practical interest for the last two decades [1-7] . Early 
direct extensions of the Chapman-Enskog method to inelastic systems have 
been refined recently so that the transport coefficients for the Boltzmann 
gas are now known accurately [5,8]. These refinements have clarified the role 
of the time dependent reference state, the local homogeneous cooling state 
(HCS), and its effects on the qualitative and quantitative characterization of 
the conditions for the validity of a hydrodynamic description. The extension 
of these results to higher densities also has been accomplished in the context 
of the phenomenological Boltzmann-Enskog equation [6]. However, more gen- 
eral conclusions regarding transport in dense granular fluids are quite limited. 
On the other hand, for elastic collisions the analysis of transport in fluids is 
quite sophisticated using, for instance, methods of linear response theory [9]. 
Such methods provide formally exact expressions for transport properties in 
terms of time dependent fluctuations in the equilibrium reference state. These 
representations have proven to be particularly useful due to the possibility 
of their direct evaluation by molecular dynamics simulation. In this context. 
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it is natural to ask to what extent the methods of Unear response theory can 

be applied to granular fluids. 

In applications of linear response for fluids with elastic collisions, the lo- 
cal conserved densities of mass, energy, and momentum play a special role. 
They represent the slow variables as the system approaches uniformity. More 
significantly, they are the hydrodynamic modes in the long wavelength limit, 
i.e the formal eigenfunctions of the Liouville operator corresponding to the 
hydrodynamic eigenvalues. Standard linear response methods consider per- 
turbations of the homogeneous state (the Gibbs ensemble) that are linear 
combinations of the conserved densities. In this way, hydrodynamic excita- 
tions are selected and, by formal manipulation of the response functions, 
Green-Kubo expressions for the transport coefficients are derived directly. A 
quite reasonable extension of this method to granular systems, using linear 
combinations of the local conserved densities as initial perturbations, has 
been described in Ref.[10]. The resulting Green-Kubo expressions are simi- 
lar to those for systems with elastic collisions, but they do not agree at low 
density with those from the Chapman-Enskog method noted above. This dif- 
ference has been made more precise by reformulating the Chapman-Enskog 
results in a Green-Kubo representation [11]. The resulting time correlation 
functions are not simply those for the fluxes of the conserved densities, as it 
is the case for elastic collisions. The reason is that the hydrodynamic modes 
for granular systems are no longer simply related to the conserved densities, 
even in the long wavelength limit. Instead, the hydrodynamic modes and con- 
served densities are conjugate sets in the sense that the former generate the 
hydrodynamic excitation while the latter measure the response. An objective 
here is to clarify this subtle difference in the context of the linearized Boltz- 
mann equation for inelastic hard particles, and to discuss the consequences 
for more general applications. 

Clearly, an important first issue is the unambiguous definition of hydrody- 
namic modes. This is accomplished by considering the eigenvalue problem for 
the operator generating the dynamics for the Boltzmann equation, linearized 
around the reference homogeneous cooling state (see Section 3 below). The 
hydrodynamic part of its spectrum is identified first in the long wavelength 
limit. Hydrodynamics refers to systematic approximations to the macroscopic 
balance equations for the local relevant densities (here our choice is the de- 
viations of number density, temperature, and flow velocity from their val- 
ues in the homogeneous state) . Consequently, we first derive these equations 
and find they are closed without approximation in the long wavelength limit. 
These asymptotic equations determine d+2 hydrodynamic eigenvalues, where 
d is the dimensionality of the system. For normal fluids these are all zero, 
while for granular systems there are three distinct eigenvalues, one of which 
is d-fold degenerate. The corresponding eigenvalues for the linearized Boltz- 
mann equation then identify its asymptotic hydrodynamic spectrum. The 
hydrodynamic spectrum more generally is defined to be the eigenvalues that 
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are continuously connected to these asymptotic values in the long wavelength 
limit. The expression "hydrodynamic mode" means the hydrodynamic cigcn- 
function with an exponential time dependence, in the appropriate time units, 
characterized by the corresponding hydrodynamic eigenvalue. A primary new 
contribution here is the calculation of these hydrodynamic cigcnfunctions in 
the long wavelength limit, showing their difference from the local conserved 
densities. 

With the hydrodynamic eigenfunctions identified, the hydrodynamic part 
of the solution to the linearized Boltzmann equation can be extracted. This 
is the dominant contribution for large times and long wavelengths. More 
specifically, this hydrodynamic part is calculated to first order in the gra- 
dients of the hydrodynamic fields and the result is used to calculate the 
average fluxes in the macroscopic balance equations at finite wavelengths. 
This provides the Navicr-Stokcs hydrodynamics for granular gases, and al- 
lows identification of the transport coefficients. The results agree with those 
from the Chapman-Enskog method [5,8] and the Green-Kubo representation 
of Ref. [11]. However, their derivation here exposes more clearly the hydro- 
dynamic eigenfunctions and their role as the origin of key differences from 
gases with elastic collisions. 

In the next Section, the nonlinear Boltzmann equation for a granular gas 
is defined, and the exact balance equations for the density, temperature, and 
flow velocity are derived from it. The homogeneous state of these equations 
corresponds to global cooling, and the corresponding homogeneous cooling 
solution (HCS) to the nonlinear Boltzmann equation is described. This so- 
lution represents the homogeneous state about which small spatial pertur- 
bations are considered. The corresponding linear Boltzmann equation and 
linear balance equations for these small perturbations are then described. 
Next, the eigenvalue problem for the generator of the time dependence in the 
linearized Boltzmann equation is defined in Section 4, and the hydrodynamic 
eigenfunctions and eigenvalues are calculated in the long wavelength limit. 
These results are used in Section 5 to derive the Navier-Stokes hydrodynamic 
equations, with expressions for the transport coefficients in terms of the lin- 
earized Boltzmann collision operator. Their equivalence with previous results 
from the Chapman-Enskog method is demonstrated, and these expressions 
arc also given a Green-Kubo representation. Finally, the analysis is recapit- 
ulated, and implications for a more general treatment of transport in rapid 
granular fluids beyond the low density context are discussed. 
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2 Nonlinecir Boltzmann Equation and the 
Homogeneous Cooling State 

The Boltzmann equation for a system of inelastic smooth hard disks {d = 2) 
or spheres {d = 3) of mass m and diameter a has the form [3,12] 

|-+v.vV(r,r,<) = J[/|/], (1) 



where J[f\f] is the nonlinear inelastic Boltzmann collision operator, 

J[f\f]= I dvT(v,Vi)/(r,v,t)/(r,vi,t). (2) 

Here ^(v, Vi ) is the inelastic binary collision operator 

f (v, vi) = a'^-i Jd&0{g-&)g-& {a-H-' - 1) , (3) 

with g = V — Vi , O the Heaviside step function, <t a unit vector joining the 
center of the two particles at contact, a the coefficient of normal restitution, 
and an operator replacing all the velocities v and vi appearing to its 
right by their precollisional values v* and , given by 

6^ V = V = V — (<T • g) cr, b„ vi = Vi = vi + {cr ■ g) cr. (4) 

For arbitrary functions g(v,vi) and /i(v,vi), it is 

j j dvi5i(v,vi)T(v,vi)/i(v,vi) = j j dvih{v,vi)T{v,vi)g{v,vi), 

(5) 

where 

T(v,vi) =a'^-iy" d&0{g-&)g-&{b^-l). (6) 

The operator bo- changes the velocities v and Vi to its right into their post- 
coUisional values v' and v'^ , 

, / 1 + , 1 + a 

= V = V — (<T • g) a, b„vi = vi = vi H — (<t • g) a. (7) 

By using Eq. (5), the equalities 



/ 



dv I mv I J[/|/] = I I . (8) 

-¥C[/], 



mv 
1 



are easily obtained. In the above expressions, p = nT is the pressure (Boltz- 
mann's constant has been set equal to unity as usual), n being the number 
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density and T the granular temperature, and C is a a nonlinear functional of 
the distribution function, 

C[/] = ^^~^r{^)^f~' 1'^^ I '^^i/(r>v,i)/(r,vi,iy. (9) 

The hydrodynamic fields are defined in the standard way in terms of 
velocity moments, i.e. 

n{v,t)= j dvf{v,v,t), (10) 

n{r,t)u{r,t)= j dvvf{r,v,t), (11) 

^n{r,t)T{r,t) = J dv^V'f{r,v,t), (12) 

where V = V — u is the velocity of the particle relative to the local velocity 
flow u. Their balance equations follow directly from Eq. (8), 

^ + V.(nu) = 0, (13) 

— +u- Vu + (mn)-^V • P = 0, (14) 

^+u- VT+^(P: Vu + V-q)+TC = 0. (15) 
at nd 

The functional giving the pressure tensor P and the heat flux q are: 

P(r,i) =m J dvVVf{r,v,t) =p{r,t)\ + J dv D(V)/(r, v, t), (16) 

q(r,t) = ^J dvV^Vf{r,v,t) = j dv S(V)/(r, v, i). 



(17) 



Here I is the unit tensor, and 



D(V)=m(vV-ll/2|), S(V) = - V. (18) 

The Boltzmann equation (1) admits a special solution /hcs describing 
the so-called homogeneous cooling state (HCS) and having the scaling form 
[3] 

/hcs(v, t) = nHVo'^{t)xHCsic), c = , (19) 

VO{t) 

where 

,o(i)=[?Z£££Wl'^' (20) 
m 
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is the "thermal velocity" and Xhcs(c) is an isotropic function of c. As in- 
dicated in Eq. (19), all the time dependence of /hcs occurs through the 
homogeneous temperature Tncsit), which obeys the equation 

^Ili^ + THcs{t%HCs{t) = 0, (21) 

with 

(1 — a^)7r^CT''~^n//Wo(i) f , f , , ^ / m 

CHcs{t) = 2r{'^+^)d J J '^*^i^^C'^(*^)^^<=''5(^i)l*^~^il • 

(22) 

Substitution of the scaling form, Eq. (19), into the Boltzmann equation leads 
to 

1 d 

2^^'-'^dv ~ J[fHCs\fHCs]- (23) 

The explicit time dependence of the HCS can be eliminated by introducing 
a dimensionless time scale s by 

where £ = [nH(T^~^) ^ is proportional to the mean free path of the particles 
in the HCS. In terms of the s variable, Eq. (21) becomes 

§-Thcs{s) = -CoThcs{s), (25) 

with 

Co = — TTT- , (26) 
vo{t) 

that does not depend on time. Equation (25) is the differential form of the 
Haff's homogeneous cooling law [13]. Similarly, Eq. (23) takes the form 

1 d 

■^Coq^ ■ {cxhcs) = J'[xhcs\xhcs], (27) 

J'[Xhcs\xhcs] = j dcif{c,ci)xHCs{c)xHCs{ci), (28) 
^(c, ci) = j d&O [(c - ci) • <t] (c - ci) ■ (T(a-2/,^i _ (29) 

The operator h~'^ in the last equation is defined again by Eqs. (4), but sub- 
stituting (v, Vi) by (c, Ci). The form of the distribution function xhcs given 
by the solution of Eq. (27) has been extensively studied. A summary of the 
results is given in Ref. [14]. 
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3 Linearized Boltzmann Equation 

Let us now consider a situation where the inelastic gas is very close to the 
HCS, so that we can write 

/(r,v,t) = .fecs(v,i)+<5/(r,v,t), |(5/(r, v, i)| « ,fecs(v, t). (30) 

Substitution of Eq. (30) into the Boltzmann equation (1 ), keeping only terms 
up to first order in Sf yields 

^6f{r, V, i) + V • V(5/(r, v, t) = K5f{r, v, t). (31) 
The linearized Boltzmann collision operator is defined by 

KSf{r.^,.t) = J[fHcs\S.f] + J[Sf\fHcs]- (32) 

Equation (31) can be written in dimensionless form, using the scaled velocity 
c defined in Eq. (19), the reduced time s introduced in Eq. (24), and the 
dimensionless length scale 

l=f (33) 
Then, the linearized inelastic Boltzmann equation becomes 

d d 

^ + c • ^ (5x = ASx, (34) 

where 

5x{l,c,s)=nH\^it)Sfir,v,t) (35) 

and 

Adx = J'[XHCs\Sx] + J'MXHCs] - 1^ • (c^x) • (36) 

In this representation, there is not any explicit dependence on time due to 
the reference state, and all the stationary methods known from the case of 
elastic collisions [9] can now be applied to the analysis of Eq. ( 34). 

We define relative deviations of the hydrodynamic fields from their values 
in the HCS by 

= = / dcSx{l,c,s), (37) 

riH J 

u;(l,s) ^ = / dccSx{l,c,s), (38) 

vo{t) J 

dii, s) ^ = jdc{^-^-l^ dx(l, c, s), (39) 

where Sy{r,t) = y{r,t) — yHCs{t) denotes the deviation of the local value 
of the macroscopic variable y from its value in the HCS. Taking velocity 



8 J. Javier Brey et al. 

moments in the Boltzmann equation (34), the hnearized balance equations 
for the hydrodynamic fields are obtained, 

^p(l,s) + |.a;(l,s) = 0, (40) 

§~s^^^' '^^Iwi'""^^' '^^Iwi'^^^' «)+'5Co(l, s)-Co m, s) + p(l, s)] = 0. (42) 
The reduced linear pressure tensor fl and heat flux are defined as 

n = ^ [^(1, s) + p(l, s)]i + J dcA{c)Sx{l, c, s), (43) 

0(1, s) = j rfcS(c)5x(l,c,s), (44) 
where A and cf) are the dimensionless forms of D and S, 



c^. s ^ o d + 2 



A(c)=cc--I, E{c) = i^c^ - j c, (45) 

The quantities H and cj) are related with the linearization of the fluxes given 
in Eq. (16) and (17) by 

Finally, the term 5Co(l) s) is the dimensionless linear deviation of the cooling 
term from its value in the HCS, 

^^0^^'^) = /'^^ / dc,XHCs{c)Sx{l,c,s)\c-c,f. (47) 

Of course, Eqs. (40)-(42) can be derived also by linearizing the nonlinear 
balance equations (13)-(15). They are the starting point for deriving hydro- 
dynamic equations, since they describe the evolution of the hydrodynamic 
fields in the linear approximation. This requires closing the equations, trans- 
forming the formal expressions of the pressure tensor, the heat flux, and the 
cooling rate into expressions in terms of the own hydrodynamic fields. These 
are the so-called constitutive relations. 
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4 Eigenvalue Problem 

In this Section wc explore solutions to the linearized Boltzmann equation 
and characterize the hydrodynamic modes. This is done through a formal 
analysis of the eigenvalue problem for the associated linear operator. Since 
the equation is linear and the collision operator A does not change the space 
variable, it is sufBcient to consider a single Fourier mode, i.e. we take 

5x(l,c,s)=e'''-'^x(k,c,s). (48) 

The linearized Boltzmann equation, Eq. (34), then becomes 

^JX=iA-ik-c) Sx, (49) 

whose formal solution is 

^X(k, c, s) = e^f^-*-) ^x(k, c, 0) . (50) 

The above expression shows most clearly that the possible excitations for the 

gas under these conditions of small initial perturbations are determined by 
the spectrum of the linear operator A — ik ■ c. This suggests a study of the 
eigenvalue problem 

{A-ik- c) c) = A, (k) e.(k, c). (51) 

The problem is posed in a Hilbert space of functions of c with a scalar 
product given by 

<g\h>= j dcxHcsic)9*{c)h{c), (52) 

where g* denotes the complex conjugate of g. Of particular interest here are 
the eigenfunctions and eigenvalues X^(c) associated with linear hydrody- 
namics. These can be identified as follows. First, at k = they correspond 
to those eigenvalues of A alone that coincide with the zeroth order in the 
gradients eigenvalues of the linearized balance equations (40)-(42), which are 
easily seen to be 0, Co/2, and —Co/2, i.e. 

^^^(c) = A^^^(c), A^=0,±|, (53) 

where £,'^{c) stands for ^^(k = O.c). More generally, at finite k the hydrody- 
namic modes are defined to be the solutions of Eq. (51) that are continuously 
connected to those of Eq. (53) as fc — > 0. The eigenvalue problem (53) is 
solved below, so the existence of the hydrodynamic modes is assured. The 
formal solution (50) then can be decomposed in the form 

<5x(k,c,s) =^S^(k,s)^^(k,c) + ^x"(k,c,s). (54) 

/3 
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The first term arises from the hydrodynamic modes while the second one 
represents all the other "microscopic" excitations. The hydrodynamic coeffi- 
cients a^(k, s) can be expressed as 



S^(k,s) = e^^^C') < ^^(k,c)|(5x(k,c,0) >=< ^"^(k, c)|<5x(k, c, s) > . (55) 



Since yl — ik-c is not self-adjoint, the functions C/j(k, c) differ from C^(k, c) but 
together they must form a biorthogonal set. It will appear below that these 
coefficients are closely related to the hydrodynamic fields of the linearized 
balance equations. 

The expected context for a hydrodynamic description is now clear. For 
given initial conditions both hydrodynamic and microscopic excitations con- 
tribute to (54). However, for long wavelengths (small k) it is expected that 
the microscopic excitations decay on time scales short compared to those for 
the hydrodynamic modes, leaving only the latter on long time scales. The rel- 
evant separation of time scales is expected to occur after a few collisions, so 
there is a large and interesting time scale for the hydrodynamic description. 
Preliminary analysis based on kinetic models for the Boltzmann equation [15- 
17] and also on low order matrix representations of the linearized Boltzmann 
collision operator [18], suggest that the separation of modes assumed here 
does occur, and hydrodynamics does dominate for wavelengths and times 
long compared with the mean free path and the mean free time, respectively. 
Furthermore, it appears from these models that the cooling rate remains 
smaller than the collision frequency even for strong inelasticity, so that inclu- 
sion of the non-conserved energy among the hydrodynamic variables seems 
justified. 

Hydrodynamic modes for fe = 

We consider now the eigenproblem associated with the linearized homoge- 
neous Boltzmann operator. 



Of course, finding all the solutions of this equation is an impossible task, even 
in the elastic limit. Nevertheless, it is quite easy to obtain some particular 
solutions, which will turn out to be the relevant ones in the hydrodynamic 
regime. We begin by constructing the function 



where C = c — w. Although F has the functional form of a "local" HCS 
distribution, the parameters p, u), and 7 do not need to have any particular 
meaning here, playing just the role of auxiliary parameters. From Eq. (27) it 
is obtained that F verifies 



(56) 



F(c,p,w,7) = pxHCsilC) 



(57) 



|Co^-(CF)=7'^+V'[F|J^]. 



(58) 
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Then, taking derivatives with respect to p, w, and 7 in the above equation 
and, afterwards, the Umit p = 7= l ,u> = 0, leads to the set of equations 

^V'l(c) = -iCo^3(c), ^V2(C) = ^CoV'2(c), ylV'3(c) = -iCoV'3(c), (59) 

where 

^i(c)^ (^|^)^ = Xi/cs(c), (60) 

(dF\ d 
^j^-XHCs(c)d = -^-[cxHCs(c)]. (62) 

Here the zero subindex indicates that the derivatives are evaluated at p = 
7 = l,u> = 0. From Eqs. (59) it is easily seen that the set of functions 

{^0{C)} = {Vl(c) - V3(C),^2(C),^3(C)} (63) 

are solutions of Eq. (56) corresponding to the eigenvalues 

Ai=0, A2 = |, A3 = -|, (64) 

respectively. The eigenvalue A2 is, therefore, d-fold degenerate. This confirms 
(53) above, indicating points in the spectrum at fc = corresponding to 
the linear balance equations. Hence, these are the hydrodynamic modes at 
longest wavelengths, C^(c). 

The eigenfunctions ^^(c) can be interpreted as follows. Consider a lo- 
cal form of the HCS distribution, /j^^^, whereby the hydrodynamics fields 
are replaced by their actual space and time dependent values in a general 
nonequilibrium state, 

/llcs(r,v,<) = n{v,t)v'i^-\t)xHCs{C), (65) 

where 



m 



(66) 



This is the analogous to the local equilibrium distribution for elastic colli- 
sions. It represents a state whose averages values of 1, v , and v'^ are those 
of the nonequilibrium state. This seems a reasonable choice to describe an 
experimentally prepared state. The linearization of such an state around the 
HCS gives 

/Hcs(r, V, t) ^ /ffcs(v, t) + nHVo\t)Sx^^\c, 1, s), (67) 
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5x'^'\c, 1, s) = [^^{c) + e^'(c)] p{l, s) + t2{c) ■ s) + (c)^?(l, s). (68) 

Thus, the eigenfunctions ^^(c) describe the hnear deviation of the local HCS, 
appearing then as the natural perturbations for a macroscopic preparation 
of the initial nonequilibrium state. This will be discussed further below. 

It is easy to show that the functions {^^(c)} are linearly independent so 
that they define a d + 2 dimensional subspace of the Hilbert space. Although 
the functions ip 13(c) arc orthogonal with the definition of scalar product in 
Eq. (52), the eigenfunctions C^(c) are not, as a consequence of the operator A 
being non-Hermitian. This leads, in principle, to consider the left eigenvalue 
problem and, therefore, to study the adjoint operator yl+ defined by 

< g\Ah >*=< h\A+g >, (69) 

for arbitrary g{c) and h{c). Using the property given in Eq. (5) it is found 

^"^5(c) = y dciXHCs{c)xHCs{ci)T{c,ci) [xhcs(c)5(c) + Xhcs(ci)5(ci)] 

+|xHCs(c)c • ^ [xhcs(c)5(c)] . (70) 

The eigenfunctions of yl+ would provide the left eigenfunctions of A. We have 
not been able to find all the eigenfunctions corresponding to the eigenvalues 
in Eq. (64). Nevertheless, for our purposes here, and also for many other 

applications, it will suffice to identify a set of functions |^^(c),/? = 1,2,3| 
verifying the biorthonormality condition 

<€"^|^^, >=5;3,/3., (71) 

/?, /3' = 1, 2, 3. A convenient choice is given by 

eT(c)=X//C5(c), 12(c) =cxffcs(c), ^3 = {j + l)xHCs{c). (72) 

The above functions span a dual subspace of that spanned by the right eigen- 
functions {^^(c)}. The fact that they are linear combinations of 1, c, and 
is closely related with the existence and approach to a hydrodynamic de- 
scription, as it will be shown in the next Section. Let us stress that while 
(c) and (c) arc eigenfunctions of A'^ corresponding to the eigenvalues 
and Co/2, respectively, the function ^3(0) is not an eigenfunction of The 
important point is that for any linear combination 

3 

ff(c) = ^a^^^(c), (73) 

the coefficients are given by 

«/3 =< t^lg >= I rfcxHCs(c)f^(c)5(c). (74) 
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It is instructive to consider what happens in the elastic Hmit a — > 1. Then 
XHCsic) becomes the MaxweUian Xmb(c) = tt""^/^ exp(— c^) and the right 
eigenfunctions are given by 

{ [2c' -{d+ 1)] Xmb(c), 2cxmb(c), {2c' - d)xMB(c)} . (75) 

Of course, these functions are not orthogonal and differ from those usually 
employed in the linear analysis of the elastic Boltzmann equation [19]. This 
is an important point, since when solving the eigenvalue problem associated 
with the inelastic homogeneous Boltzmann equation by means of a pertur- 
bation calculus around the elastic limit, the first step in the perturbation 
scheme is to remove the degeneracy of the elastic zeroth eigenvalue by exactly 
solving the perturbed eigenvalue problem within the subspace spanned by the 
degenerate eigenfunctions. The relevant conclusion is that the functions ^^(c) 
appear in a natural way even when considering the small inelasticity limit. 

5 Navier- Stokes and Green-Kubo expressions 

In the last section the hydrodynamic modes were defined and calculated at 
asymptotically long wavelengths. It remains to identify the hydrodynamic 
coefficients in Eq. (54) as well as the hydrodynamic modes at finite k. Con- 
sider first the coefficients defined by Eq. (55 ) and evaluate the biorthogonal 
function ^^(k, c) to leading order in k, so that 

a'^ik,s)-^<^^ic)\Sx{k,c,s)> . (76) 

Use of the expressions given in Eq. (72) shows that, in coordinate represen- 
tation, these are just the hydrodynamic fields of the linear balance equations 

{a^(I, s)} ^ |p(I, s)Mi, s), ^9{l, s) + p{l, s)| . (77) 

This completes the connection between the hydrodynamic modes and the 

dynamics determined from the linear balance equations. They are the same 
for long wavelengths and small perturbations, as expected. 

5.1 Navier-Stokes Approximation 

To obtain the hydrodynamic eigenvalues (k) to Navier-Stokes (fc^) order, 
the eigenvalue problem (51) can be solved treating ic • k as a small pertur- 
bation of the operator A [20]. This method implicitly restricts k to be small 
relative to the deviation of the restitution coefficient a from unity, the elastic 
limit. A more uniform description, without restrictions on k relative to a, is 
obtained by calculating the fluxes in the balance equations as functions of 
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a^(k, s), to first order in k. The resulting Navier-Stokes hydrodynamic equa- 
tions can then be used to determine (k). This is tlie approacli used liere. 
Tlie first step is accomplislied by using tlie identity 

valid for arbitrary operators A and B, in the formal solution (50) to get 

5x(k, c, s) = e"^5x(k, c, 0) - / ds' e^'^ik ■ ce-^'^'^-^ '^^^xCk, c, s). (79) 

Jo 

Next, substitute the representation (54) to identify the dependence on a|g(k, s), 

(5x(k,c,s) = e"'^(5x(k, c, 0) - / ds' e^'^ik ■ ce-^'^^-*-'^) 

Jo 

^S;^(k,5)C^(k,c)+^x'"(k,c,s) 
^ e«^<5x(k, c, 0) - V [ rfs'e'*'^ce-*'^^^(c) ■ikS^(k,s) 
= e«^5x(k,c,0)-^ [rd5'e^'(^-^^)c€^(c) 



• ika^(k,s). 



(80) 

In the second transformation, it has been assumed that k is sufficiently 
small so that only leading order terms need to be retained, and that s is 
sufficiently large to neglect the microscopic excitations relative to the hy- 
drodynamic modes contributions. Moreover, in the last equality use has been 
made of the fact that is an eigenfunction of A. In coordinate representation 
this result becomes more explicitly 



where 
and 



fi 

Sx(}, c, s) ^ e'^Sxih c, 0) - ^ F0{c, s) ■ s), 

/3=1 



{1^0(1,3)} = {p{l,s),u:{l,s),e{l,s)} 

Fi(c,s)= / dsV^a(c)c + 2F3(c,s), 

Jo 



(81) 

(82) 

(83) 

(84) 
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F3(c,s) 



(85) 



The pressure tensor and the heat flux up to first order in the gradients of 
the fields are now obtained by substituting Eq. (81) into Eqs. (43) and (44). 
The term involving the initial condition (5^(1, c, 0) docs not contribute to these 
quantities since, in the long time limit, it is given by a linear combination 
of the functions x^(c), that are all orthogonal to both A{c)xhcs{c) and 
IJ{c)xHCs{c)- Taking into account the symmetry properties of the system, 
the resulting expressions can be written in the form 



niiii,s) = -[e{i,s) + pii,s)]5ij 



dh 



LOj{l,s) 



dl 



-LOiil, s) 



2d_ 
ddi 



■u}(l, s)Sij 



0(1,S) = -K^0(1,S) 



(86) 



(87) 



Equation (86) is the expected Navier-Stokes expression for the pressure ten- 
sor, involving the shear viscosity coefficient 77, but Eq. (87) contains, besides 
the usual Fourier law characterized by the heat conductivity k, an additional 
contribution proportional to the density gradient and with an associated 

transport coefficient /i. This latter term has no analogue in elastic gases. 
The expressions of the (time-dependent) transport coefficients are 



1 f 
= + / ^c/iij(c)F2,y(c,s) 

k{s) = ^J dcS{c)-F3{c,s), 



fj,= 



dcl!{c) ■ Fi(c, s). 



(88) 



(89) 



(90) 



As already mentioned, the fimctions Aij{c)xHCsic) and Si{c)xHCsic) are 
orthogonal to the subspace spanned by the right eigenfunctions ^/3(c). There- 
fore, the decay of the "correlation functions" in Eqs. (88) (90) lies outside the 
spectrum of the slowest modes, and the s-timc integral in the expression of 
the functions converges. Consequently, it is possible to set s — > 00 for times 
scales long compared with the microscopic times. In this way, the reduced 
transport coefficients become timc-indcpendcnt. Returning to the original 
variables, the above expressions for the hydrodynamic fluxes are equivalent 
to 



nnTHCsit) 
d 



1 



6T{r,t) 6n{r,t) 



Tncsit) nn 



-V 



_d_ 



Sui 



-(V-,5u)4 



(91) 
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q(r, t) = -K V5T(r, t) - /x V^n(r, t) , (92) 
with the transport coefficients given by 

■q = nHmlva{t)fi, (93) 

K = nHivQ(t)k, (94) 

M = ^T^M- (95) 



5.2 Green-Kubo Relations 



It is possible to write Eqs. (88)-(90) in the form of Green-Kubo relations. 
Define an operator A by 



j dcg*{c)Ah{c) = j dch*{c)Ag{c), 



(96) 



for arbitrary h and g. Comparison with Eq. (69) yields 

Mc) = Xffcs(c)^+ [9{c)XHCs{c)] 

= J dcixi/cs(ci)r(c, ci)[<?(c) + 5(ci)] + • ^g{c). (97) 

Also define "correlation functions" in the HCS by 



<gh>= I dcxHCs{c)g{c)h{c). 



Then, the expressions for the transport coefficients can be rewritten as 



1 f°° 

K=- ds < E(s)-^3 > e"'^''/^ 
dJo 



jl = 2k + 



dJo 



ds < I!{s) ■ *i >, 



with the definitions 



^2,»j(c) 



1 



#1 = C - 2*3(c), 

dlnxHCs{c) 



*3(c) = 



d + c 



' da ' 

91nxHCs(c) 



dc 



(98) 

(99) 

(100) 
(101) 

(102) 

(103) 

(104) 
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The time dependence of the dynamical variables is given by 

g{c,s)=e'^g{c). (105) 

Equations (99)-(101) are the Green-Kubo formulae for a dilute granular gas 
described by the Boltzmann equation. They express the transport coefficients 

as integrals of time-correlations functions in the HCS. Although they are not 
identical to those derived in Ref. [11] by a different method, it is easy to show 
they are equivalent, since they differ in terms giving vanishing contributions. 

Let us also compare them with the expressions derived from the nonlinear 
Boltzmann equation by the Chapmann-Enskog method [5,8]. Consider the 
reduced shear viscosity coefficient. It is given by Eq. (88) with F2Aj defined 
in Eq. (84). In the limit of large s, we have seen that the integral in the latter 
expression can be extended up to infinity. Performing then the integral, it is 
obtained 

or, equivalently, 

where is the solution of the integral equation 

(^.l-|)c,,(c) = -6,i(c)c,-, (108) 

being orthogonal to the subspace spanned by the lowest order right cigcn- 
functions of A. Equations (107) and (108) are equivalent to those derived by 
the Chapman-Enskog procedure in Refs. [5,8], where the explicit expressions 
of the transport coefficients in the first Sonine approximation is obtained. 
A similar analysis leading to the same conclusion can be carried out for the 
other transport coefficients. The only conceptually relevant difference is in 
the reference state. In the present linearized analysis the transport coeffi- 
cients are defined in the HCS, whereas in the Chapmann-Enskog method 
they are defined in a local HCS, characterized by the parameters n{r,t), 
u(r,f) and T{r,t). As a consequence, these are the fields appearing in the 
scaling in Eqs. (93)-(95). 

It is also instructive to consider the elastic limit a — » 1 of the Green- 
Kubo expressions derived here. In that case, we already know that the HCS 
distribution becomes the Maxwellian Xmb, and Co — * 0. Then, the functions 
defined in Eqs. (102)-(104) become 

*i(c) ^ -2S{c), ^2,y(c) ^ 2ciCj, *3(c) -> U - ^] c. (109) 
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Therefore, the transport coefficient jl, given by Eq. (101) vanishes in the 
elastic hmit, as expected, and 



2nmivo{t) 
(P + d-2 

2nmivo{t) 
d? + d-2 



y2 ds < Aij{s)ciCj > 
V / ds < Aij{s)A,j > 



(110) 



^^n,= '^ r ds<E{s).U-'^]c> 



' ds < Uis) ■E>. (Ill) 



d 



f 

Jo 



The last equalities in the above equations follow from the symmetry of the 
time correlation functions. For elastic systems, the temperature is constant 
to lowest order in the gradients of the hydrodynamic fields and, therefore, in 
the integrands appearing in the above expressions the definition of the scale 
s, given by Eq. (24), reduces to s = Vot/£, i.e. it is simply proportional to 
the actual time. Then, Eqs. (110) and ( 111) are the low density limit of the 
usual Green-Kubo expressions for the transport coefficients in terms of the 
autocorrelation functions of the microscopic fluxes of the relevant densities 
[9]. 



6 Discussion 

The primary results of the above analysis have been the identification of 
hydrodynamic modes in the spectrum of the linearized Boltzmann operator, 
their leading order calculation in Section 4, and their use to derive the Navier- 
Stokes approximations to the heat and momentum fluxes. Of particular in- 
terest are the long wavelength forms of the hydrodynamic eigenfunctions 
given in Eq. (63), which are not simply linear combinations of the conserved 
densities l,v,v^. It is clear from (54) that perturbations of the homogeneous 
state due to ^^(c) generate the hydrodynamic excitations, while the response 
to these excitations is measured by the fields a^{s). These fields are averages 
of the biorthogonal set of functions , which are linear combinations of the 
conserved densities 1, v, i;^, as required by the macroscopic balance equations. 
As anticipated in the Introduction, this "conjugate" relationship between the 

sets and {C/?} is a key difference between normal and granular gases. 

For elastic collisions, the two sets of flelds are the same. This has impor- 
tant consequences for the derivation of transport properties as it has been 
illustrated here. For example, the reduction leading to the Navier-Stokes ap- 
proximation in Eq. (80), depends critically on the fact that the functions 
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are eigenfunctions of A. Any other choice for the initial perturbation would 

lead to a different time dependence, mixing microscopic and hydrodynamic 
excitations, and precluding direct identification of transport properties. 

The identification of the hydrodynamic modes discussed here is important 
for the application of linear response methods [21] to granular systems. In 
order to illustrate it, let us consider again the formal solution to the linear 
Boltzmann equation given in (50) and choose for the initial condition the 
linearization of a local HCS as given by Eq. (67), 

^x(k,c, s) = e^c^-'"--) J2 < ^pimo) > e^(c) 

13 

= e^(^-^''-)^a0(k,s = O)C^(c). (112) 



This initial state clearly represents a hydrodynamic perturbation. The re- 
sponse a/3(k, s) is then given by 

a/siKs) =J^Ri30'{k,s)ai3'{k,s = 0), Ri30'{k,s) =< ^^le^^^-*'^)^^, > . 
0' 

(113) 

The response functions i?/3/3'(k, s) are time correlation functions comprised of 
the pair of the conjugate densities. It can be shown that application of stan- 
dard linear response methods to these response functions leads to the correct 
Green-Kubo relations with the proper convergence properties for the time 
integrals. Conversely, the same procedure applied to correlation functions 
comprised of the {^0} alone does not lead directly to well-posed Green-Kubo 
expressions. Let us also mention that linear perturbations of the kind consid- 
ered in Eq. (112) have been used [22]in direct Monte Carlo simulations of the 
Boltzmann equation [23], to verify the validity of the linear hydrodynamic 
description and the explicit expressions of the transport coefficients. 

Although the analysis here has been limited to the low density linearized 
Boltzmann equation, these considerations point the way for a proper applica- 
tion of linear response methods in the more general context of nonequilibrium 
statistical mechanics. 
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